Dynamic transition and Shapiro step melting in a frustrated Josephson-junction array 
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We consider a two-dimensional fully frustrated Josephson-junction array driven by combined di- 
rect and alternating currents. Interplay between the mode locking phenomenon, manifested by 
giant Shapiro steps in the current-voltage characteristics, and the dynamic phase transition is in- 
vestigated at finite temperatures. Melting of Shapiro steps due to thermal fluctuations is shown to 
be accompanied by the dynamic phase transition, the universality class of which is also discussed. 
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Driven systems with many degrees of freedom arc ubiq- 
uitous in nature and exhibit a rich variety of dynamic 
behavior as the driving or the temperature is varied. 
Here the usual equilibrium concepts are not applica- 
ble, making the study of nonequilibrium transitions in 
such driven systems very important and highly nontriv- 
ial. A prototype system is provided by the Josephson- 
junction array under driving currents ^^^^ which is 
also closely related to, e.g., moving vortex lattices and 
sliding charge density wavesi^ In equilibrium the fully 
frustrated Josephson-junction array (FFJJA), as an ex- 
perimental realization of the fully frustrated XY model 
with the U(l) x Z2 symmetry, has been a source of contro- 
versy as to whether a Z2 transition and a U(l) transition 
occur at the same temperatureAiSiiiiiS In the presence 
of direct driving currents, the transition temperatures in 
general reduce^ approaching zero as the current is in- 
creased toward the critical current; beyond it, the Z2 
transition occurs at a substantially lower temperature 
than the U(l) transition. 4 On the other hand, when the 
driving current has only the ac component, the FFJJA 
has been shown to exhibit a dynamic Z2 transition, which 
belongs to the same universality class as the equilibrium 
Z2 transition^ 

The remaining case is the FFJJA driven by combined 
direct and alternating currents, which is well known to 
display quantized voltage plateaus, called giant Shapiro 
steps, in the current-voltage characteristics at zero tem- 
perature^* This is a manifestation of mode locking in 
response to the driving current, associated with topolog- 
ical invariance of the system. 14 However, the issue of the 
FFJJA driven by combined direct and alternating cur- 
rents has not been addressed from the viewpoint of the 
dynamic phase transition, which requires analysis of the 
system at finite temperatures. 

This work investigates how such a mode locking phe- 
nomenon changes into the dynamic transition, as the 
temperature is raised from zero. To describe the dy- 
namic transition associated with topological invariance, 
we introduce the dynamic order parameter and investi- 
gate its behavior with the driving amplitude. The zero- 
temperature behavior of the dynamic order parameter is 
discussed in the framework of the topological invariance 



of Shapiro steps. As the temperature is raised from zero, 
the Shapiro steps are shown to melt, which is accompa- 
nied by the dynamic phase transition. From the detailed 
behavior of the dynamic order parameter at finite tem- 
peratures, we construct the dynamic phase diagram on 
the plane of the temperature and the direct current, and 
examine nature of the transition. 

To begin with, we consider the equations of motion 
for phase angles {</> z } of the superconducting order pa- 
rameters in grains, which form an L x L square lattice. 
Within the resistively-shunted-junction model under the 
fluctuating twist boundary conditions, 15 they read: 



d<t>jj 
dt 



sin(0jj - Tij ■ A) + ijij 



0, (1) 



where the primed summation runs over the nearest neigh- 
bors of grain i and the thermal noise current Tjij satisfies 
{Vij{t)Vki(t')) = 2T5(t—t')(8ikSji — 8u5jk) at temperature 

T. We have used the abbreviations 4>ij = 4>i ~ 4>j ~ Aij 
and = 1*3 — Vj with = (xj, yt) denoting the position 
of grain i. Note that for nearest neighboring grains 
becomes a unit vector with the lattice constant set equal 
to unity. We have also written the energy and the time in 
units of Hi c /2e and H/2eRi c , respectively, with the criti- 
cal current i c and the shunt resistance R of a single junc- 
tion. The dynamics of the twist variable A = (A X1 A^) 
is governed by the equations 

= ^2 sin (^j ~ + ^A* ~ r dc ~ he sin fit 



dAy 

dt 



(2) 



where X^i-A denotes the summation over all nearest 
neighboring pairs in the fi (= x, y) direction, r]& sat- 
isfies (vA,{t)VA^{t')) = {2T/L 2 )S(t~t')5^,, and the 
combined direct and alternating currents Idc + lac sin fit 
are injected in the x direction. In the Landau gauge, the 
bond angle Ay, given by the line integral of the vector 
potential, vanishes for Vj = r j + x and takes the value 
TTXi for r , = r, + y. 
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At zero temperature the FFJ JA driven by combined di- 
rect and alternating currents displays integer, fractional, 
and subharmonic giant Shapiro stepspi&A^ The system 
evolves periodically in time with the topological invari- 
ant period, characteristic of each voltage step. Time evo- 
lution of the phase configuration shows that the system 
visits periodically various accessible states consisting of 
the ground states, the transient states, and other dynam- 
ically accessible states of the purely dynamic origin>i& To 
characterize such dynamic behavior associated with the 
Z 2 symmetry in the FFJJA, we consider the chirality 



C(R, t) = sgn | sin [<fe (*) ~ r « ' A W] | 
and the staggered magnetization 



(3) 



(4) 



R 



where ^ p denotes the directional plaquette summation 
of links around the dual lattice site R = rj + (l/2)(x+y). 

Figure [2 exhibits the zero-temperature time evolution 
of the staggered magnetization on various voltage steps, 
where V denotes the time-averaged voltage in units of 
Ltth/2e. The period of m(t) in Fig.QJis given by r, 2r, 
and 4r on steps V = 1, 1/2, and 1/4, respectively. Here 
we have considered the FFJJA of size up to L = 24, 
and integrated numerically the equations of motion in 
Eqs. ||TJ and J5J via the modified Euler method, using 
time steps of size At = 0.01. The values of At have 
been varied, only to give insignificant difference. Typi- 
cally, data have been averaged over 1000 driving periods 
of the (alternating) current, after the initial 500 periods 
discarded; this has turned out sufficient for reaching ap- 
propriate stationarity. 

At high temperatures thermal fluctuations are so 
strong that the influence of the driving current and the 
lattice potential can be neglected, leading to random fluc- 
tuations of m(t). We thus expect a phase transition be- 
tween the dynamically ordered phase and the disordered 
one as the temperature is varied. To describe such a dy- 
namic phase transition, we introduce the dynamic order 
parameter, defined to be the staggered magnetization av- 
eraged over n periods of the (alternating) current: 



Inn 



t +nr 



m(t) dt 



(5) 



with t ee 2ir/n. On each Shapiro step, we consider the 
topologically invariant period^ to choose the proper in- 
teger n. Specifically, on the l/2g step we choose n = q, 
avoiding the null order parameter. In general, the dy- 
namic order parameter in Eq. J2J depends on to and 
henceforth the maximum value is defined to be the dy- 
namic order parameter: Q = max to Q(ta). 

Figure |21 shows the ensemble average (Q) of the dy- 
namic order parameter as a function of the temperature 
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FIG. 1: Time evolution of the staggered magnetization m(t) 
at zero temperature. The amplitude and the frequency of 
the ac component are I ac = 1.0 and Q./2-k — 0.1, respectively. 
The dc component is given by Idc = 0.05, 0.177, 0.40, and 0.75 
from the top to the bottom, corresponding to the average 
voltage V = 0,1/4,1/2 and 1. The period of m(t) is r (= 
2ir/Q = 10), 2r, and 4r on the steps V = 1,1/2 and 1/4, 
respectively. For clarity, data for V = and 1/4 have been 
shifted upward, while those for V = 1 shifted downward. 
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FIG. 2: Dynamic order parameter as a function of the temper- 
ature in the system of size L = 24 for I dc = 0.04(D), 0.4(O), 
and 0.75(A), corresponding to V = 0, 1/2, and 1, respectively. 
Lines are merely guides to eyes. Inset: Time averaged volt- 
age V versus the dc component Idc at temperature T = 0, 0.1, 
and 0.24. As T is raised, the voltage plateaus tend to melt. 



T at various values of the dc component Idc- Manifested 
is the presence of a dynamic phase transition, separat- 
ing the dynamically ordered phase at low temperatures 
and the disordered phase at high temperatures. In all 
the three cases in Fig. |2 corresponding to the voltage 
steps V — 0,1/2, and 1, the dynamic order parameter 
increases gradually as the temperature is lowered from 
high temperatures. Subsequently, it grows rapidly, even- 
tually saturating to the zero-temperature value. Further, 
it is also observed that cooling and heating curves for the 
dynamic order parameter exhibit appreciable hysteresis 
at low temperatures (not shown here), which arises from 
the asymmetry of the lattice potential induced by the 
dc component of the driving current and the resulting 
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FIG. 3: Dynamic phase diagram on the T—Id c plane. The 
phase boundaries are determined by the crossing points of 
Binder's cumulants for size L = 8, 16, and 24. The areas en- 
closed by the boundaries specify dynamically ordered phases 
characterized by Q ^ 0. Lines are merely guides to eyes. 



anisotropy, similarly to the dc driven case^ Such behav- 
ior does not appear in the system driven by alternating 
currents only (without the dc component). 

The Shapiro step originates from the interplay between 
the period of the external driving current and the pe- 
riod emerging from the intrinsic dynamics of the system. 
Such mode locking behavior present in the low temper- 
ature regime is expected to smear out as the thermal 
fluctuations become strong. Indeed the Shapiro steps 
are observed to melt as the temperature is raised, disap- 
pearing eventually at sufficiently high temperatures (see 
the inset of Fig. [2J . To probe the relation between this 
melting behavior and the dynamic phase transition, we 
first estimate the transition temperature accurately via 
Binder's cumulanljii 



U T =1- 



(Q A ) 

3<Q 2 r 



(6) 



which takes the value 2/3 at T = 0, vanishes in the high- 
temperature limit, and becomes size- independent at the 
transition temperature. The unique crossing point of Ul, 
independent of the system size L, for given I e i c thus yields 
the corresponding data point for the phase boundary on 
the T—Idc plane in Fig. |3 The areas enclosed by the 
resulting phase boundaries specify the dynamically or- 
dered phase, characterized by Q ^ 0, for V = 0, 1/2, 
and 1 from the bottom to the top. The transition tem- 
perature initially decreases monotonically to zero as the 
driving amplitude Idc is increased in the region of V = 0. 
Further increase of Idc drives the system into another or- 
dered region, characterized by V = 1/2, where the tran- 
sition temperature first grows with the driving current, 
then reduces to zero, and so on. 

For comparison, phase boundaries are also estimated 
independently from melting of the Shapiro steps. Here we 
define the width of a step as follows: For a given point on 
a voltage step, we measure the slope of the straight line 
connecting the point to the middle point of the step and 
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FIG. 4: (a) Width w of the voltage plateau versus the tem- 
perature T on the step V = 0(D, +), 1/2(0, x), and 1(A, V), 
where the two symbols for each step represent the data from 
two different criteria of the slope: 0.005 and 0.01, respectively, 
(b) Dynamic phase diagram on the T—Idc plane, determined 
by the width of voltage plateaus for size L = 16. The areas 
enclosed by the phase boundaries specify dynamically ordered 
phases with w 7^ 0. Lines are merely guides to eyes. Inset: 
Detailed diagram displaying the ordered phase corresponding 
to the subharmonic voltage step V = 1/4, again determined 
by its width. 



consider the point to be located on the step if the slope 
is less than 0.005. Namely, the step width is determined 
by the boundary point, for which the slope is given by 
0.005 [the use of other values hardly changes the result- 
ing phase boundaries, as shown in Fig. Ufa)]. The step 
width obtained in this manner diminishes with the tem- 
perature and eventually vanishes to zero [see Fig. Ufa)], 
from which each phase boundary in Fig. 0fb) has been 
estimated. The good agreement between Figs.|3and0Jb) 
manifests the correspondence between the nonzero value 
of the dynamic order parameter and the finite width of 
the Shapiro step, thus indicating that the dynamic phase 
transition accompanies melting of Shapiro steps. It is 
also expected that there exist additional dynamically or- 
dered phases, associated with the series of subharmonic 
Shapiro steps^ they presumably form Arnold tongue 
structure^ in the complete phase diagram. We have thus 
examined the subharmonic step V = 1/4, and indeed 
found the corresponding ordered phase, shown in the in- 
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FIG. 5: Scaling plot of the dynamic order parameter for size 
L = 8(D), 16(0), and 24(A) at dc component I dc = 0.04. 
The fitting has been performed with the critical exponents 
v = 0.82 and f3/v = 0.11. 

set of Fig.^Jb). Unfortunately, other phases, correspond- 
ing to higher-order subharmonic steps, occupy even tinier 
regions, making it very difficult to locate the boundaries. 

We finally study the nature of the transition by means 
of finite-size scaling for the dynamic order parameter: 

(Q) = L-^fdT-TjL 1 /"). (7) 

In Fig. H we plot (Q)L^_ U versus (T - T C )L X I V at I dc = 
0.04 (corresponding to V = 0), which yields reasonable 
scaling collapse with the critical exponents v = 0.82 and 



f3/v = 0.11. Since these values of the critical exponents 
agree well with those for the equilibrium Z 2 transition 
in the FFXY model^ill we conclude that the dynamic 
phase transition in the system here belongs to the same 
universality class as the equilibrium Z2 transition. Simi- 
lar conclusion was also reached in the FFJ JA under weak 
staggered oscillating fields or uniform alternating cur- 
rentsi^iii This suggests that the universality class of the 
transition does not change if V remains to be zero (i.e., in 
the same Arnold tongue). On the other hand, at larger 
values of 1& C corresponding to V ^ 0, the scaling plot 
does not collapse well and there still exists ambiguity on 
the nature of the transition; this appears similar to the 
dc driven case, where the transition nature changes as 
the driving current exceeds the critical value^^ 

In summary, we have examined the 2D FFJJA driven 
by combined direct and alternating currents, as a proto- 
type example of driven systems exhibiting rich dynamic 
behaviors. The relation between melting of the voltage 
plateaus (Shapiro steps) and the dynamic phase transi- 
tion has been explored. It has been observed that melt- 
ing of the Shapiro steps is accompanied by the dynamic 
phase transition, revealing the correspondence between 
the two phenomena. We have also examined the nature 
of the transition, to find, on the zero voltage step, the 
same universality class as the equilibrium Z2 transition. 
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